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Abstract 



We study the relaxation of interacting single-molecule magnets (SMMs) in both spatially ordered 
and disordered systems. The tunneling window is assumed to be, as in Fes, much narrower than 
the dipolar field spread. We show that relaxation in disordered systems differs qualitatively from 
relaxation in fully occupied cubic and Fes lattices. We also study how line shapes that develop in 
"hole-digging" experiments evolve with time t in these fully occupied lattices. We show (1) that 
the dipolar field h scales as t v in these hole line shapes and show (2) how p varies with lattice 
structure. Line shapes are not, in general, Lorentzian. More specifically, in the lower portion of 
the hole, they behave as (| h | / i? 5 ) 1 if h is outside the tunnel window. This is in agreement 
with experiment and with our own Monte Carlo results. 

PACS numbers: 75.45. +j, 75.50.Xx 
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I. INTRODUCTION 



Magnetic relaxation in crystals of single-molecule magnets (SMM's), such as Fes, has 
become a subject of great interest.— At low temperature, T, each SMM behaves ap- 
proximately as a single spin S. Magnetic relaxation at k^T < 0.1U/S, where U is a mag- 
netocrystalline anisotropy barrier, is temperature independent, and is duly attributed to 
quantum tunneling under the barrier. Hyperfme interactions with nuclear spins as well as 
dipole-dipole interactions among all SMM electronic spins give rise to a variety of phenom- 
ena that are not yet fully understood. Hyperfme interactions enable spins to tunnel even 
when the ensuing Zeeman energy change 2Eh is much larger than the tunnel splitting energy 
A, provided | | is smaller than some e w .^ For Feg, for instance, A ~ 10~ 4 mK, e w ~ 10 
niK, and the rms value of the Zeeman energy 5eh is approximately 400 mK. We shall restrict 
ourselves to systems with e w <C 5eh- In these systems, the relaxation of the magnetization 
goes on long after time T -1 , where V is the spin tunneling rate for spins within the tunneling 
energy window (that is, spins for which | Sh \< s w ). Tunneling spins give rise to changing 
dipolar fields, which in turn bring new spins into the tunneling energy window, thus keeping 
a magnetic relaxation process from extinction. In this paper, we focus our attention on 
effects that stem from this process. 

We consider here experiments of the sort that were reported by Wernsdorfer et al. in 
Ref. {§?]. In them, a crystalline sample of SMM's is first quenched from k^T ~ U/S to 
ksT < 0.1U/S in either (a) a weak applied magnetic field of a few mT, and later observed 
after the field is switched off at time t = or (b) a zero field, and later observed after a weak 
field is applied at t — 0. We shall refer to the former as a FC (for field cooled) experiment 
and to the latter as a ZFC (for zero field cooled) experiment. Both types of experiments 
can be lumped into one defining fa = 1 — m/m , where mo is the initial magnetization in a 
FC experiment, and m = m/m s , where m s is the final steady state magnetization in a ZFC 
experiment. Wernsdorfer et al. observed the m oc t 1 ! 2 over roughly two time decades. The 
time evolution of "holes" that, under suitable conditions, develop in a magnetization density 
function, have also been reported.^ The existing theory at the tim o 5 i 8 i 9 predicted a universal 
\/t short time relaxation from a fully polarized system, but said little about relaxation from 
or into weakly polarized states.— Other theories^ take hyperfme interactions into account 
but disregard dipole-dipole interactions. They therefore apply if 5eh < e w which is not 
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within the scope of this paper. 

We have developed a theor y 12 ' 14 that gives the time evolution both of fh and of the 
holes' line shapes in weakly polarized systems of interacting SMMs, such as Fe§, in which 
e w <C 5Eh- There are three clearly discernible time regimes. For Ft < 1, fh oc Yt. In the 
second time stage, when 1 < Yt, up to some time before fh ~ 1, fh oc t p , at least for all 
fully occupied cubic lattices. Moreover, the theory gives a simple relation that specifies 
how p varies with lattice structure. In FC experiments, fh ~ 1 in the third time stage, 
that is, m m 0. The third time stage is more interesting in ZFC experiments. Then m(t) 
settles down temporarily to a quasi stationary value m s , which the theory predicts, if either 
/cbT ^> e w or if the heat exchange rate with the lattice is much smaller than T; on the other 
hand, if k^T 3> e w is not fulfilled and heat exchange rate with the lattice is not much smaller 
than T, then the relaxation of the magnetization shifts into a thermally driven approach to 
equilibrium, skipping the quasi stationary state. A quite different treatment of relaxation 
from weakly polarized states that gives 1 — m/rriQ oc \/t, independently of the spins' spatial 
distribution, is given by Tupitsyn, Stamp, and Prokof'ev (TSP) in Ref. [l^, criticized in 
Ref. [3], and defended in Ref. Q]. TSP's treatment of relaxation from weakly polarized 
states is rather unrelated to the earlier theory^ of Prokofe'v and Stamp for a \ft relaxation 
from fully polarized states, about which we have nothing to say here. According to TSP, 
^/t relaxation in weakly polarized systems holds as long as 1 < Yt and fh -C 1, that is, 
roughly, over the second time stage. This is also the time domain where our theory gives 
fh ~ (e w /5eh)(Yt) p for fully occupied lattices. This time stage, sometimes referred to as 
"short times", can in fact be arbitrarily large, for arbitrarily small e w , since m <C 1 only 
implies Yt <C (fe^/e^) 1 ^. Note also that \ft relaxation from fully polarized systems ends 
sooner, when t ~ 5eh/e w ^ This may explain why Wernsdorfer— et al. were able to observe 
fh ~ \fi in Fe§ (p ~ 0.5 for Fe§) under weak applied magnetic fields for up to some 10 3 s 
while Ohmii et al. observed \/t relaxation from fully polarized Fes for only up to 10 2 s. 

The two approaches, ours and TSP's, are quite different. The underlying assumptions 
follows. In Ref. |l5j, the main result follows from a claim that is made on the function 
f(h,t) = p^(h,t) — p^(h,t), where p^(h,t) and Pi(h,t) are the number densities of up- and 
down-spins, respectively, with a dipolar field h acting on them at time t. In its latest form,— 
the claim is that the scaling form f(h, t) = f(h/t 2 ) "was found to be valid in our MC simu- 
lations for different lattice types..." when 1 < Yt and m <C 1. Since, the magnetization m 
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and f(h, t) are clearly related by m{t) = — J dhf(h, t), a \/t relaxation follows immediately. 
As we have shown recently— the above scaling form holds approximately for SC and Feg 
lattices (as defined in Ref. la]), but not in general. It fails, for example, for FCC, BCC, 
and diamond lattices.— Our theory also gives the time evolution of f(h,t), but follows from 
a more fundamental assumption: that the dipolar field on any one given site changes by 
some random amount Ah, whenever a spin flips somewhere else for the first time, and that 
Ah follows a Lorentzian distribution (for more details see Sects. Ill Al and IIII Al and Ref. 

). The integro-differential equations for the evolution of f(h, t) and of the magnetization 
that obtain in our theory, follow from this assumption. 

Unfortunately, as far as we know, only experiments on a crystalline Fes structure have 
thus far been performed. However, MC simulations have been performed for various fully 
occupied cubic lattices, which have given values of p that agree with our predictions . 13 ' 
This paper's first aim is to extract from our theory how the magnetization is supposed to 
relax in Fe§, compare this with experiment over the time span where published experimental 
data exists,-^ and make predictions for later times. It is also our purpose to predict, and 
check with MC simulations, how m relaxes with t for other spatial distributions, namely, 
under full spatial disorder, thus providing another test for our theory. 

The holes observed in the experiments described above& are also of interest. They cor- 
respond to "wells" that develop in the function f(h,t) = pi(h,t) — p^(h,t) (defined above). 
From the relation m(t) = — J dhf(h, t), the time evolution of m follows, but f(h, t) provides 
additional information about the magnetic evolution of the system that m(t) does not.^ For 
short times, that is, for Ft < 1, the hole's width is equal to s w Q This was first surmised 
by Wernsdorfer et to propose a number, approximately 10 mK, for the tunneling energy 
window, e w . However, we know of no published data for the holes' line shape evolution well 
into the intermediate time range, that is, for 1 <C Ft. Our second aim is to fill this gap. 
To this end, we work out from our theory the time evolution of holes' line shapes in this 
time stage in fully occupied cubic systems and Feg crystals, and check the results obtained 
against our MC results. We also obtain new results of a more general nature for the holes' 
line shape. Before we state our results, we specify the model. 
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A. The model 



All spins are on a lattice, they point along the easy anisotropy axis, and interact as 
magnetic dipoles. We consider both fully and partially occupied lattices. Let the magnetic 
field at site i produced by spin Sj at site j be given, in the usual notation, by = 
hd{a/rij) 3 [l — 3(zij/rij) 2 }, where is the distance between the i and j sites, a is the 
distance between nearest neighbor sites, h d = (fi /A7i)gfi B S/a 3 , and Si = ±S for all i. 
Furthermore, let the magnetic field hi at site % be given by hi = ^ • hy, where ^ ■ is over 
all occupied lattice sites. The tunnel window size and tunneling rate F are defined next. At 
very low temperature, that is, if k B T < 0.1U/S, a spin can flip only if the field h acting on 
it satisfies | h \< h w M The flipping rate is F if upon tunneling the energy decreases, but if 
the energy increases by AE, then the rate is rexp(— AE/k B T), following detailed balance. 
(Even though | AE \< e w , and ususally k B T 3> e w , exp(—AE/k B T) is not quite equal to 
1, and this makes a difference after a sufficiently long time, as is shown below.) We also 
simulate relaxation processes in which the energy E is assumed to remain constant (such 
as if no spin-lattice relaxation takes place). Then, we assume a value of T is such that E 
remains approximately constant. No tunneling window restriction applies for spin flips if 



We let p(h, t) be the probability density function (PDF) that any one given spin have field 
h at time t, let po(h) be the same distribution for a completely random spin configuration, 
and let 



For a Gaussian field distribution, a is equal to the dipole field rms value Sh for a random 
spin configuration (see Table I), but this is not so in general. Values of a that follow from 
MC simulations for cubic and Fe§ lattices with randomly oriented spins are given in Table 



Randomly oriented spins on a cubic lattice give a Lorentzian field distribution of ho half 
width at half maximum if n lM Values we will be using for ho are given in Table I. From 
here on, unless otherwise stated, all magnetic fields and energies are given in terms hd and 
gfiBh d S, respectively. 



k B T > U/S. 




(1) 



I. Finally, let h = (8tt 2 /3 5/2 )h d n, where n is the number of dipoles per unit cubic cell.— 
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TABLE I: Quantities ho, cr 7 and yj (/i 2 )o are given for randomly oriented spins on the lattices 
specified. h = n8vr 2 /3 5/2 , o = [\/2vfp(0)] _1 , where p(0) stands for the PDF at h = 0, and y/(h 2 } 



is the rms spatial average of the dipolar field. For all lattices except for Fes, ho, a, and W (/i 2 )o are 
given in terms of h&. All lattices are fully occupied, except for the "dilute SC", which stands for a 
SC lattice with h <C 1 sites occuppied. Finally, a = 87r 2 /3 5 / 2 . 



LATTICE 



o 



a 



SC 
BCC 
FCC 
Fe 8 6 
dilute SC 



a 

2a 

4o 

47(1) mT 

ah 



3.655 
3.864 
8.303 
46(1) mT 



24 J 



3.83(2) 
4.03(2) 
8.44(2) 
31(1) mT 



We assume an easy anisotropy axis as given in Refs. |23j|2j|2 



B. Plan and main results 



Section |H] is devoted to the relaxation of the magnetization. The equations from our 
theory^ which we use to calculate the time evolution of the magnetization are restated in 
Sect. Ill Al The results that follow from them for Fes crystals are shown to agree rather 
well with experiment and with our own MC results in Sect. Ill Bl The evolution we predict 
for m(t) as well as our MC results cover a time span that is 2 orders of magnitude longer 
than the the experimental one te {te ~ 20 minutes). Let t w be the end time for the regime 
where m oc t p . We showed in in Ref. (l4J that t w w T~ 1 (a/h w ) 1 / p , from which we obtain 
t w ~ IOte. For t w < t, the evolution of m(t) is shown to depend sensitively on T if good 
thermal contact with a heat reservoir is assumed. If on the other hand we assume constant 
energy processes (i.e., no spin-lattice relaxation), then m levels off, if only temporarily, to a 
stationary value m s when r w < t. The value of m s we obtain from theory is unrelated to the 
equilibrium value of m, which only obtains much later. This stage, when m m s , sets in 
after most spins in the system have tunneled at least once after the magnetic field is applied. 
Simulations bear this out. We also obtain, from theory as well as from MC simulations, m(t) 
for spatially disordered systems. More specifically, we make a random selection of a fraction 
n of L x L x L SC lattice sites and place spins on them. For n < 0.1, we assume full disorder. 
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Then, theory predicts a magnetic relaxation that bears no resemblance to a \ft rule, not 
even to the t p rule that we obtain for fully occupied SC lattices. Instead, for n < 0.1, 

a 

approximately, where r m ~ 10 6 r _1 and q ~ —0.105. Results from our MC simulations are 
in fair agreement with this. In Sect. 11111 we report results for holes' line shapes in fully 
occupied crystal lattices. The main results for line shapes, which are derived in Sect. 1111 Al 
follow. When 1 « Tt « {(y/h w f/ p and h w <| h + H |< a, 

f(h,t)/f(h,0)^\r ] \ 1 P ~ 1 (3) 

if | 77 |<C 1, where 

= k + H (A) 

H is an applied field, and p is given by 

p^ 1 sin up = \p2^Kajh^. (5) 

In Sect. IIII Al we also derive relation for holes line shapes that hold over longer time spans: 
1 ^ Ft < [o le w ) 1 l p . In Sect. IIII Bl we apply these results to published^ experimental data 
for Fes and to MC data for Fes as well as to fully occupied FCC lattices. Finally, concluding 
remarks appear in Sect. IIVI 



II. RELAXATION OF THE MAGNETIZATION 



A. Theory 

We first describe a stochastic model which helps to understand the physics of the problem 
as well as the statistical assumption we have made in order to solve it. Consider two tracks, 
both filled with particles. Let there be one particle on the "up-track" for each up spin on 
a lattice, and one particle in the "down-track" for each down spin. Let all particles on 
each track be ordered according to the value of the magnetic field H + h acting on each 
spin. To mimic tunneling, random select a particle within the tunnel window, that is, a 
particle fulfilling —h w <H + h<h w , whether on the up- or down-track, and move it to the 
"point" H + h on the opposite track. In order to mimic the effect such a spin flip has on the 
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dipolar fields on other spins, draw a random value of Ah for each particle from a Lorentzian 
distribution of half width ho/N, where iV is the total number of spins, and let h — > h + Ah 
if the particle that just shifted track has done so for the first time. This latter proviso is 
related to the fact that no effect on the di pola r field follows when the same spin flips twice 
(for a more detailed explanation, see Ref. [.14]). Repeat this whole process at every tic of a 
clock. Clearly, the whole process stops when all particles haved jumped track at least once. 

We can get a feeling for the relevance of lattice structure or spatial spin distribution 
from the following simple consideration. Np(0)2h w is the number of particles in the tunnel 
window, which is the number of times the clock tics in time which, in turn, would 
be the relaxation time for the equilibration of the number of up and down particles in the 
tunneling window if there were no "field shifting". Now, the median field shift, or diffusion 
length, in time T _1 is 2p(0)hoh w , that is, a fraction p(0)ho of the tunnel window's width. 
Thus, p(0)ho is a measure of the amount by which the unbalance between up and down 
particles is restored in the tunnel window in time T -1 . Therefore, the relaxation rate clearly 
depends on p(0)ho, which in turn depends on lattice structure. This shows why the latter 
is relevant to the relaxation of the magnetization. 

The simple statistical assumption above has enabled us to derive^* the equations we need 
for the calculation of the time evolution of the magnetization. These equations, which we 
reproduce in a compact form immediately below, give the magnetization m(t) at time t, 
and n(t), the fractional number of spins that flip at least once in time t. We first recall an 
important ingredient of the theory for ZFC experiments^ the energy per spin at the time 
when the system is quenched, which we refer to as the "annealing energy", is — e a . Let X\ = 
mgfi B Sa{h 2 ) /(e a h w H) and X\ = 2(a/h w )(l—m/m ), for ZFC and FC experiments (defined 
in the Introduction), respectively, and X2 = na/h w for both FC and ZFC experiments. The 
desired equation follows, 



a\ = 4, h\ = 2. In order to obtain x\, Eq. © must first be solved for j = 2, letting a% — 1 
and &2 = 1, in order to then use x 2 (t) in Eq. (J7|), and thus enable substitution of u into Eq. 
© for j = 1. 





where 





S 



The theory applies if h w <C cr and the energy is constant, that is, if no energy tranfer 
to the phonon bath takes place. This is also approximately so if kT 3> e w . If the constant 
energy condition is not met in a ZFC experiment, a linear in time magnetization relaxation 



that is thermally driven takes over before m(t) —>■ r?v— The theory also gives 



2He a /(gSfi B (h 2 ) ). 



(8) 



Note that the definition of Xi(t), together with Eqs. (jHJ), (JZJ), and (jHI) imply that the time 
variation of m(t)/m s in a ZFC experiment is the same as 1 — m/mo in a FC experiment. 
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FIG. 1: to and n (the fractional number of spins that have flipped in time t at least once) versus 
time. Full lines are for theory and symbols stand for MC results for 65536 spins that are initially 
up or down at random, with probabilities 0.6 and 0.4, respectively, on a FCC lattice. ♦ (o) stand 
for to for h w = 0.1 (0.01); □ stand for n for h w = 0.01. The dash-dotted lines stand for slopes 
p = 0.73, and p = 0.5, which follow for FCC lattices from Eq. © and from Ref. [j ], respectively. 
Data points taken from Ref. 3] are shown as (). Our data points for Tt < 1024 come from 
averages over 1400 MC runs, and for 100 runs for t > 1200. 



A few remarks about Eqs. © and (J7J) are in order. Clearly, Xi(t) only depends on 
two parameters: h /cr and cr/h w . The latter only comes into play at the later portion of 
the time evolution, when n(t) > 2a 2 / '(vt/iq) , which is when m(t) starts leveling off. Now, 
0.4 < cr/h < a/tt/2 for all cubic lattices, whether fully occupied or not. It follows that 
leveling off of m(t) is triggered some time (see below) while n(t) < 1. This has to do with 
the fact that a spin flip contributes to hole widening in f(h,t) the first time it takes place 
after the field is switched on. When a spin flips a second time, it only returns to its initial 
state, thus it cancelling the effect os the first flip. Since m is proportional to the width 



of the hole in f(h,t) when Ft > 1 it follows that m(t) becomes constant when n ~ 1. 
The time when m(t) levels off is illustrated in Fig. for fully occupied FCC lattices. For 
instance, for h w = 0.01, theory predicts m(t) to start leveling off when n(t) ~ 0.11 since 
2<t 2 /(7t/zq) ~ 0.11 then. Thus, a significant portion of the evolution, up to n ~ 1, only 
depends on h /a. Information about the lattice or degree of spatial disorder only comes 
into the equations through this number. Other numbers, such as h w , e a , and H only come 
into the proportionality factor between x\ and either m or 1 — m/rriQ. The temperature does 
not enter anywhere into the equations. 

We have shown analytically in Ref. ^| that m oc t p and n oc t p fulfill Eq. (JHJ) for j = 1 
and j = 2 in the 1 <C Tt <C (<y/h w ) 1 / p time range if p is given by Eq. (j5J). Numerical 
solutions of Eq. (jUJ) shows that m oc t p ensues in the wider 1 < Tt < time span 

for all fully occupied cubic lattices.— We show below that this is also so for a fully occupied 
Fe 8 lattice. 

Spatially disordered systems behave differently. Consider a very small fraction (n <C 1) of 
sites of a SC lattice to be occuppied by randomly oriented spins. Then, a Lorentzian dipolar 
magnetic field distribution^ of half-width ho ensues. It follows then, from the definition of 
cr and of h 0l that <j/h = a/ it/2. The number p = follows then from Eq. (jSj). Recall, 
however that theory only implies that this ensues only when Tt ^> 1. For earlier times, more 
specifically, for all Tt > 1, we find that the numerical solution from Eq. (jUJ) is well fitted 
by Eq. (J2J) if n < 0.1. Numerical solutions of Eq. © are plotted in Fig. for n = 1 and 
n = 0.6, both for h w = 0.02 and in Fig. |2b for (1) n = 0.1 and h w = 0.02, and for (2) 
n = 0.03 and h w = 0.006. Note that the same solution obtains for the latter two cases. We 
come back to these figures in Sect. Ill Bl 



B. Comparison with experiments and simulations 

We first make use of Eqs. © and (|7j) to obtain m(t) for Fes- Some numbers must first 
be fed into Eqs. (Hp and ([7|). For h w , we usq2£ 0.8 mT, as given in Ref. We use T = 0.04 
s _1 , (see Refs. |l4| and j^). With the numbers given in Table I for a and h , we obtain 
Xi(t) and X2(t) numerically from Eqs. (JB1) and ((Zj). Finally, the value of — e a , the annealing 
energy— is needed in order to obtain m from X\. Not knowing e a , we treat it as a fitting 
parameter. We find e a ~ 36 mK fits best the experimental data points from Ref. [6J), which 
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FIG. 2: (a) Magnetization versus Ft. Lines are for theory and symbols are for data from MC 
simulations of L x L x L SC lattices with a fraction h of their sites occupied. For h = 0.6 the 
dashed line is for theory, 0, and ♦ are for MC data for L = 32 and L = 16, respectively. For n = 1 
the full line is for theory, A, and ▲ are for MC data for L = 32 and L = 16, respectively. All 
symbols are from averages of 800 MC runs while Ft < 10 3 and 40 MC runs when Ft > 10 3 . (b) 
Same as in (a) but for n = 0.1 and h w = 0.02, and for h = 0.03 and h w = 0.006. The full line is 
for theory for both values of n. For h = 0.03 o, and • stand for MC data for L = 64 and L = 32, 
respectively. For n = 0.1, t>, and ► stand for MC data for L = 32 and L = 16, respectively. As in 
(a), all symbols are from averages of 800 MC runs while Tt < 10 3 and 40 MC runs when Ft > 10 . 

are shown in Fig. [3] for a few applied fields. The energy —36 mK may be compared to the 
approxiamte value, —500 mK, of the ground state energy.— 

The MC data points shown in Fig. El follow from simulations in which the system first 
evolves at some high temperature (a few Kelvin) for a short time (less than 1 MC sweep) until 
the energy equals —36 mK. At such temperatures, Fe§ cluster spins are not forced to tunnel 
through the ground state doublet. Accordingly, all spins are allowed to flip, regardless of the 
dipolar field acting on them. We explore different scenarios after quenching. In our theory, 
we assume no energy exchange takes place between the spin system and a heat reservoir. We 
have also performed MC simulations under this assumption. This is approximately realized 
for the time range exhibited in Fig. (jSJ) by flipping only spins within the tunnel window, 
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and then with equal probabilities for upward or downward flips. If, on the other hand, 
heat exchange does take place readily, as one g 

rates that are comparable to T are found, then detailed balance should be enforced in MC 
simulations. The results of doing this lead to the plots shown in Fig. |3] for T = 40 and 300 

HlK. 



a 




t[h] 

FIG. 3: Magnetization versus time, in hours. Units for m are such that m = 1 for full polarization. 
0, □ and o are for experimental data taken from Ref. [6j for H = 3.92, 2.24, and 1.12 mT, 
respectively. Full lines are from our MC simulations for kT 3> £ w , and dashed lines are for 
theoretical predictions, that is, from Eqs. © and 0- Data points that follow from MC simulations 
for H = 3.92 mT are also shown for T = 40 (♦) and T = 300 mK (A). We assume d 6 4 4 i 21 h w = 0.8 
mT and used the values of a and ho that are given in Table I. In the simulations, the initial state 
was prepared at T = 2 K. At this temperature we allowed the simulation to proceed in time up 
to the point when the energy of the system reached —36 mK, which is 0.07 of the ground state 
energy. This is the value of e a we used in order to relate x\ and m, just above Eq. (JHJ). We treat 
one MC sweep as Tt = 1 and assume^ T = 0.04 s _1 in order to convert MC sweeps to hours. Note 
that t w ~ 3 hours, since a ~ 31 mT and h w ~ 0.8 mT. 

Results obtained from theory, in Sect. Ill Al for very disordered systems are shown in Fig. 
Eh for n = 0.1 and h w = 0.02 and for n = 0.03 and h w = 0.006. Monte Carlo data points 
are also shown for the same values of n and of h w . Data points for n — 0.1 and h w = 0.02 
fall on top of data points for n = 0.03 and h w = 0.006. This is as expected, since a oc n, 
for for full spatial disorder, implies that a/h w has the same value in both cases and theory 
predicts independence from any other parameter. The fitting function from Eq. (j2J) falls 
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right on top of the curve for theory in Fig. |2b, and cannot therefore be shown separately. 

Finally, we consider size effects. For a 16 x 16 x 16 lattice and n = 0.1, for instance, 
approximately 410 spins make up the system. Of these, only approximately a fraction 
2p(0)h w are within the tunnel window. That is, approximately 820h w /irh spins, which only 
amounts to some 12 spins, are within the tunnel window. For this reason, we also simulated 
32 x 32 x 32 and 64 x 64 x 64 lattices for n = 0.1 and n = 0.03 respectively. Monte Carlo 
data points are shown in Fig. |21 Clearly, no significant size effects are observed. 

III. TIME EVOLUTION OF THE LINE SHAPE 

In this section we first derive some results for f(h,t) that are valid whenever m oc t p 
holds. We know from Ref. [14J and from the previous section that m oc t p holds in the time 
span 1 < Ft < (a/h w ) l l p for fully occupied SC, FCC, BCC, and Fe 8 lattices, but we now 
know it is not so for spatially random systems. The results we derive below are applied to 
fully occupied Fes and FCC lattices and are compared to results from experiment (for Fes) 
and from MC simulations. 



A. Theory 



The starting point for the derivation are the following two equations, from Ref. 



14] 



g(h,t)~ \ dT^^G(h + H,t-r) (9) 



o 



and 



where g(h, t) = f(h, 0) — f(h, t), u(t — r) = h^nit — r), and n(t — r) is the fractional number 
of spins that flip at least once in time t — r. The rationale for these two equations is given 
next, but (if | h + H |< a and Tt < (a/h w ) 1 / p ) Eqs. Q and flUD also follow from Eqs. 



(13) and (14) of Ref. [lj], respectively. Assume that, between times t and r, a fraction 
n(t — t) of all spins flip at least once and that n(t — r) ^ 1. This can later be checked 
to be fulfilled if t < t w , where r w = T-\a /h w ) l / p . Then, Eq. fljft gives the probability 
density G(h, H,t — t) that, at time t, the field is h + H at a site where the field at time r 
was 0.— To understand Eq. ©, note first that the definition of g(h,t) implies that g(h,t) 
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FIG. 4: (a) / versus rj, defined in Eq. @, for the times shown. All symbols stand for Wernsdorfer's 
experimental data 6 for Fes f° r I h \> 2h w . Thin lines are for data points for | h \< h w , to which Eq. 
(j3J) does not apply. The full line stands for Eq. (fT2|) . using h w = 0.8 mT and p = 0.58, from Eq. 
®, and a /ho = 0.66, from Table I. (b) Same as in (a) but for MC simulations of Fes- Points from 
within the tunnel window have been excluded. Times are given in seconds, after converting each 
MC sweep into 25 seconds, which is the value we assign to r . All symbols stand for averages 
over 1.6 x 10 5 MC runs of systems of 16 x 16 x 16 spins. The full line stands for Eq. (|12|) . 



must fulfill f dhg(h,t) = m(t) — m(0). Equation (jUJ) does give m{t) — m(0) = J drdm/dr, 
since J dhG(h, H,t — r) = 1. Similarly, a variation in the magnetization (dm/ dr)dr coming 
from some spin nipping between times r and r + dr, when the field h + H acting on them 
was within the tunnel window, contributes to g(h,t) with a Lorentzian curve whose width 
is h n(t — r), where n(t — t)/2 is approximately the fraction of the total number of sites 
where spins at time t point opposite to the way they did at time r.— Furthermore, the area 
under the Lorentzian must be given by (dm/ 'dr)dr. That explains Eq. (JHJ). 
We make use of n(t — r) ~ m(t — r),— and of 

m ( T ) ~ l.l^(rr) p (11) 
a 

which has been shown to hold for all fully occupied cubic and Fe 8 lattices in Ref. [3| and 
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in Sect. Ill Bl respectively, while 1 < (Tt) p < a/h w . We then divide by /(/i,0), which, for 
h <C a is approximately given by itlq/ ' \phxo .— Finally, the change of variable x = rjt brings, 
if 1 < Tt < (a/h w ) l / p and ft,^ <C | h + H | <C a, Eq. © into 

^ —2 I n 3^. ( 12 ) 



/(/i,0) tt 7 m? 2 + (1 -x) 2 ^ 

where a ~ 0.8(<r//io) 2 - Both Eqs. (J3J) and (JT2J) show that the field /i scales as t p in holes' 
line shapes. 

Finally, to obtain Eq. © from Eq. (JT2J), note first that /(0, t)//(0, 0) -> as t -> oo 
in Eq. (|12jl . since the integral therein equals 7r/ sinp7r if 77 = 0. 30 Then, breaking up the 
integration interval into two pieces, (1) from to (ar/ 2 ) 1//2p , and (2) from (ar] 2 ) 1 ^ 2p to 1, and 
expanding x p /(arf + x 2p ) in powers of e and 1/e in the first and second integration intervals, 
respectively, where e = x 2p /ar] 2 , gives | 77 j 1 /^ 1 for the leading term, which is desired result, 
that is, Eq. ©. 



B. Comparison with experiments and simulations 

In this section we test our results, that is, the validity of Eqs. (jSJ) and (|12p. against 
experiments^ and against our MC simulations. 

We first apply Eqs. © and (f*H2|) to Fe 8 . From Table I, a /ho = 0.66 follows. Substitution 
of this number into Eq. © gives p = 0.58. Knowing the value of p enables us to plot the 
data points for Fes shown in Fig. |3k- Unfortunately, data for holes in Fes have only been 
published for t < 40 s, that is, for Ft < 1.6, a time which falls short of the validity range for 
Eqs. (J3J) and (JTJj). Still, one can appreciate in Fig. how the data points seem to approach 
the theory curve for f(h,t)/f(h,0) as t increases up to Tt < 1.6. 

In order to see how this would go for longer times, we have used our model to simulate 
an experiment on Fes- The results are shown in Fig. Hfc- We have let one MC sweep equal 
Tt = 1, which, by the argument given above, implies t ~ 25 s for Fe§. The agreement with 
theory is remarkable. This is better appreciated in the log-log plot shown, with the same 
data, in Figs. and b. On the other hand, rescaling these plots, using p — 1/2 gives rise to 
some data point scatter, but not sufficiently large to convincingly rule out p = 1/2. This is 
not too surprising, given the small difference between p = 1/2 and the value p = 0.58 that 
is given by Eq. (jSJ). Still, one might have hoped that these data would have been sufficient 
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FIG. 5: (a)Everything is as in Fig. @k, except that a Log-log scale is used here, (b) Same as in 
(a) but for MC simulations, instead of experiments. Everything else is as in Fig. ^p. The dashed 
line stands for rj 1 ^" 1 , as predicted by Eq. (j3J) for p = 0.58, for Fes- 

to discriminate between Eq. (jSJ), where rj is raised to the 1/p—l power and the Lorentzian 
curve of Ref. Q|. Again, data for smaller values of rj would be required for this. We know 
of no other experimental results for hole digging we can make use of. As far as we know, all 
other reported experiments for SMM systems start from strongly polarized initial states.— 
Consequently, we decided to do simulations of SMM's in FCC lattices, because Eq. (J3J) 
gives then a value, 0.73, for p, which differs significantly from 1/2. We are now at liberty 
to choose the value of h w . In order to be able to obtain holes' line shapes down to rather 
small values of r], and still meet the validity criterion for Eqs. (jBJ) and (jEJ), we let h w take 
values down to 0.01. 

We show how m varies with t in Fig. [I]in a simulated experiment in which all spins are 
up and down with probabilities 0.6 and 0.4 in the initial state, and evolve thereafter with no 
applied field. For Ft < 1 (not shown), rh oc Ft. Note that m oc (Tt) p up to Tt ~ (a/hw) 1 ^, 
which for h w = 0.01, for instance, Tt ~ 9 x 10 3 , as predicted.— Note also the good agreeement 
with the value, 0.73, which Eq. © gives for p. 
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FIG. 6: (a) f(h,t)/f(h,0) versus rj, denned in Eq. (JIJ, for the shown times. Times are in MC 
sweeps. h w = 0.01. Symbols stand for averages over 1400 MC runs for 65536 spins in a FCC 
lattice. The full line is from Eq. (|12[). As in Fig. 4, initially, all 65536 spins are randomly up 
or down with probabilities 0.6 and 0.4, respectively, (b) Same as in (a) but in a Log-Log scale. 
The solid line stands for Eq. (|12|) and the dashed line is for the best fitting Lorentzian curve. The 
dash-dot line stands for rj 1 / p ~ 1 , as predicted by Eq. ©, for p = 0.73, given by Eq. © for FCC 
lattices. 

Holes' line shapes obtained from MC simulations are shown in Figs. and b. The nice 
agreement with theory is reassuring. Similar plots but using p = 0.5 give unsatisfactorily 
wide data point scatter. The data clearly follow Eq. (JHJ) for 77 <C 1 and deviate sharply 
from a Lorentzian line shape. We do not exhibit results for SC or BCC lattices, but we have 
found them to follow our predictions equally well. 



IV. CONCLUSIONS 

Results that follow from our theory for the relaxation of the magnetization of interacting 
SMMs are reported. They are in fair agreement with the experimental relaxation of the 
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magnetization observed in Fes as well as with our own MC results for Fes and for other 
lattices. Furthermore, we make some predictions for Fes that can be checked experimentally. 
Experiments following the lead of Ref. [f| would be made feasible by the application of a 
transverse field H± of approximately 0.3 T, since T, which increases as A 2 , would then 
increase by a factor of roughly 50 if H± is applied along the easier magnetization direction 
on the xy plane (see Figs. 2 and 3 of Ref. |4j). This would in effect approximately reduce 
the time scale in Fig. El from hours to minutes. Comparison of experimental results with 
MC data shown in Fig. El would be interesting. It would, for instance, show whether 
heat exchange takes place readily, as inferred in Ref. [2^ from nuclear magnetic resonance 
experiments, or not. 

We have also shown, from theory and MC simulations, that the magnetization of spatially 
disordered SMMs relaxes, not as any power of time, but approximately as given by Eq. © . 
A counterintuitive prediction that follows from our theory and from MC simulations can 
be gathered from Figs. and EJx One might have thought that dilution would lead to 
weaker dipole interactions and, consequently, to unhindered, faster relaxation. Instead, the 
opposite effect takes place for 0.1 < h < 1 after some time. 

Line shapes that develop in crystals of Fes clusters have been obtained from our theory. 
We have shown that f(h,t) is only a function of h/t p for all | h \> h w , that is, for all h 
outside the tunnel window. This is the main content of Eq. (fT2l . Furthermore, we have 
shown that data points from experiments on Fes, taken from Ref. p|, as well as results from 
MC simulations we have performed for the same system, follow this rule. Scaling also ensues 
for the data from our MC simulations of SMM's on FCC lattices for p = 0.73, as given by 
Eq. but not for the otherwise predicte d 13ll5lls p =1/2 value that is supposed to hold 
universally. 

We have also shown that f(h,t) ~| h/t p (wp)- 1 if h w /a <C i] <C 1 and h is outside the 
tunnel window. Again, this is in agreement with experimental and MC results for Fes (see 
Figs. EK and b), FCC (see Figs. Ek and b), and (not shown) SC and BCC lattices. A rough 
argument that explains why holes' line shapes are not Lorentzian follows. Note first that 
while field distributions from dilute systems of dipoles are indeed Lorentzian,— only spin 
flips that take place after time t contribute to the diffusion of a hole that was "dug" at time 
t. Since a full hole is only dug gradually in the course of time, a sum of Lorentzian functions 
of h [see Eq. (|T2J)] of various widths is expected. Not surprisingly, Lorentzian function does 
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not ensue for f(h,t) [see, Eq. (J3J)]. Here, experimental data for holes in the hundreds of 
seconds time range, over which the magnetization has already been observed experimentally, 
would be helpful. 
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